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Abstract 

There  has  been  a  recent  interest  in  modelling  the  transient  behaviour  of  proton  exchange  membrane  ( PEM )  fuel  cells.  In  the  past,  there 
have  been  several  electrochemical  models  which  predicted  the  steady-state  behaviour  of  fuel  cells  by  estimating  the  equilibrium  cell  voltage 
for  a  particular  set  of  operating  conditions.  These  operating  conditions  included  reactant  gas  concentrations  and  pressures,  and  operating 
current.  Steady-state  behaviour  is  very  common  and  in  some  cases  is  considered  as  the  normal  operating  standard.  Unsteady-stale  behaviour, 
however,  is  becoming  more  of  an  issue,  especially  for  the  transportation  applications  of  fuel  cells  where  the  operating  conditions  will  normally 
change  with  time.  For  example,  system  start-up.  system  shut-down,  and  large  changes  in  the  power  level  may  be  accompanied  by  changes  in 
the  stack  temperature,  as  well  as  changes  in  the  reactant  gas  concentrations  at  the  electrode  surface.  Therefore,  both  mass  and  heat  transfer 
transient  features  must  be  incorporated  into  an  electrochemical  model  to  form  an  overall  model  predicting  transient  responses  by  the  stack.  A 
thermal  model  for  a  Ballard  Mark  V  35-cell  5  kW  PEM  fuel  cell  stack  has  been  developed  by  performing  mass  and  energy  balances  on  the 
stack.  The  thermal  characterization  of  the  stack  included  determining  the  changes  in  the  sensible  heat  of  the  anode,  cathode,  and  water 
circulation  streams,  the  theoretical  energy  release  due  to  the  reaction,  the  electrical  energy  produced  by  the  fuel  cell,  and  the  heat  loss  from 
the  surface  of  the  stack.  This  thermal  model  was  coupled  to  a  previously-developed  electrochemical  model  linking  the  power  produced  by  the 
stack  and  the  stack  temperature  to  the  amount  and  method  of  heal  removal  from  the  stack.  This  electrochemical  model  calculates  the  power 
output  of  a  PEM  fuel  cell  stack  through  the  prediction  of  the  cell  voltage  as  a  complex  function  of  operating  current,  stack  temperature, 
hydrogen  and  oxygen  gas  flowrates  and  partial  pressures.  Initially,  a  steady-state  overall  dynamic  model  (electrochemical  model  coupled  with 
the  thermal  model )  was  developed.  This  was  then  transformed  into  a  transient  model  which  predicts  fuel  cell  performance  in  terms  of  cell 
voltage  output  and  heat  losses  as  a  function  of  time  due  to  various  changes  imposed  on  the  system. 

Keyword*:  Proton  exchange  membrane  fuel  cells;  Cell  voltage;  Heat  losses 


1.  Introduction 

The  Ballard  Mark  V  5  kW  stack  used  in  this  study  consists 
of  35  proton  exchange  membrane  (PEM)  Tuel  cells  stacked 
in  series.  The  fuel  cells  consist  of  a  DuPont  NalionIM  1 17 
proton  exchange  membrane  embedded  with  a  platinum  cat¬ 
alyst  and  with  an  active  cell  area  of 232  cm2.  The  stack  weighs 
approximately  43  kg  and  measures  approximately  38  cm  X 
21  cmX  21  cn  including  the  active  stack  with  the  internal 
humidification  section  and  cooling  cells  [  1 ) .  A  schematic  of 
the  Ballard  Mark  V  stack  [2]  is  shown  in  Fig.  I.  The  scries 
of  experimental  runs  performed  on  the  stack  were  done  with 
air  and  hydrogen  (although  the  model  has  been  developed  to 
deal  with  a  reformate  gas  feed  containing,  mainly,  hydrogen 
and  carbon  dioxide).  Extra  dry  hydrogen  from  compressed 
gas  cylinders  and  dry  purified  air  from  an  air  compressor 
enter  the  stack  separately  and  are  transported  directly  to  the 
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Fig.  1.  Schematic  of  the  Ballard  Mark  V  PEM  fuel  cell  stack  (2):  (I) 
graphite  flow  Held  plates;  (2)  hydrogen  gas  supplied  to  the  anode;  (3) 
membrane  electrode  assembly;  (4)  air  supplied  to  the  cathode,  and  (5) 
Single  fuel  ceU, 
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humid ification  section  where  the  streams  pick  up  water 
vapour  and  are  then  individually  fed  into  each  of  the  fuel 
ceils  in  their  respective  anode  and  cathode  channels  in  the 
active  fuel  cell  section.  A  circulating  liquid  v.alcr  stream 
keeps  the  stack  humidified  and  is  used  to  control  the  stack 
temperature. 


2.  Steady-state  electrochemical  model 

A  steady-state  electrochemical  model  predicting  the  volt¬ 
age  output  of  a  Ballard  Mark  V  35-cell  stack  was  previously 
developed  by  our  group  [3-5],  The  fuel  cell  voltage  was 
defined  as  the  sum  of  three  terms:  the  thermodynamic  poten¬ 
tial,  E,  the  activation  overvoltage,  rjaa,  and  the  ohmic  over¬ 
voltage,  Tjuhmii;: 

I",; cl!  f  TJnhmii: 

where 

I7api“  £|  +  £>Ts,uck+  ^Tsuh*  In  (Co2)  +  &Tslai;k  In  (0 

and 

w  -i7r=  -«(6+&raaclt+$J()  (i ) 

where  Tsmi.k  is  the  stack  temperature  (K),  i  is  the  operating 
current  (A),  cQ*  is  the  concentration  of  oxygen  at  the  cata¬ 
lyst  interface  (mol /cm3),  and  the  f  s  represent  parametric 
coefficients  based  on  experimental  data  which  can  vary  in 
value  from  stack  to  stack.  The  parametric  coefficients  &  to 
in  the  activation  overvoltage  term,  t^c„  are  based  on  the¬ 
oretical  equations  from  kinetic,  thermodynamic,  and  electro¬ 
chemistry  fundamentals.  The  parametric  coefficients  to  & 
in  the  internal  resistance  term,  /?inl,  are  purely  empirical  based 
on  temperature  and  current  experimental  data.  A  set  of  coef¬ 
ficients  for  the  Ballard  Mark  V  stack  used  in  this  study  is 
shown  in  Table  1 . 

The  effects  of  mass  transport  arc  incorporated  into  several 
of  the  three  terms  in  Eq.  ( 1 ).  The  thermodynamic  potential 
is  a  function  of  the  stack  temperature,  and  the  partial  pressures 
of  oxygen  and  hydrogen  at  the  catalyst  interface.  The  acti¬ 
vation  overvoltage  is  a  function  of  the  stack  temperature,  the 
operating  current,  and  the  oxygen  concentration  at  the  catalyst 
interface. 


experimental  cell  voltage  of  the  Ballard  Mark  V  5  kW  fuel  cell  slack. 

Fig.  2  shows  a  parity  plot  of  the  model-predicted  experi¬ 
mental  cell  voltage  versus  the  experimental  average  cell 
voltage.  It  shows  the  validation  of  the  steady-state 
electrochemical  model  using  the  coefficients  in  Table  1, 


3.  Steady-state  thermal  model 

Mass  and  energy  balances  around  the  fuel  cell  stack  were 
performed  to  calculate  the  various  energy  terms  associated 
with  fuel  cell  operation: 

4lhci>  Asciis  +  4d„,-+‘?i  llji.i  (2) 

The  energy  terms  include  the  sensible  heat  calculations  for 
each  of  the  fuel  cell  streams  ( anode,  cathode,  and  water) ,  the 
theoretical  energy  produced  by  the  fuel  cell  reaction,  the 
electrical  energy  output,  and  the  heat  loss  from  the  surface  of 
the  stack.  These  equations  are  defined  in  Fig.  3  with  symbols 
defined  in  the  List  of  symbols.  Fig,  4  shows  a  schematic  of 
the  inlet  and  outlet  streams  in  a  fuel  cell  system,  including 
the  humidification  section,  identifying  various  parameters 
including  flowrates  from  Eqs.  (3)-(9)  and  energy  terms 
from  Eqs.  ( 10)— ( 16). 


Table  I 

List  of  parametric  coefficients  for  the  Ballard  Mark  V  5  kW  stack 


Coefficient 

Estimate  for  the  activation  overvoltage  term 

Estimate  for  the  internal  resistance  term 

fi 

-0.944 

6 

3.S4X10'1 
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Fig.  3.  Thermal  analysis  equations  For  the  Ballard  Mark  V  5  kW  stack. 

The  average  heat  transfer  coefficient  for  the  stack  may  be 
estimated  using  the  average  heat  loss  from  the  surface  of  the 
fuel  cell  stack:  •,  Tain 


(  )  SIBOlC. 


(Tsucfc  —  Tmom) 


Experimental  values  of  heat  toss  were  obtained  from  the 
Ballard  Mark  V  fuel  cell  data  to  yield  an  average  ( hA )  ,lKk  of 
17W/K. 

For  comparison,  the  predicted  heat  loss  by  natural  convec¬ 
tion  can  be  estimated  using  the  Nusselt  equation  [6] : 

h  =  ja(GrPr)m  (18) 

where  a  and  m  are  numerical  parameters  which  are  found  in 
Ref.  [6] .  Note  that  the  Nusselt  equation  includes  the  Grashof, 
Gr,  and  Prandtl,  Pr,  numbers. 

Our  calculated  values  for  h  from  this  natural  convection 
correlation  range  from  about  3  to  7  W/(m2  K)  for  stack 
temperature  from  about  25  to  85  °C,  corresponding  to 
(M)sllll,k  values  from  about  1  to  2  W/K.  The  present  exper- 
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Fig.  4.  Schematic  of  the  inlet  and  outlet  streams  in  a  fuel  cell  system  iden¬ 
tifying  various  parameters  including  flowrates  and  energy  terms  defined  in 
Fig.  3. 
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Table  2 

Heat  transfer  parameters  for  the  Ballard  Mark  V  fuel  cell  stack 


Heal  transfer  parameter 

Approximate  value  ( W/K ) 

UlA), 

2 

( hA  >,■ 

10 

(M  >„. 

so 

{liA 

17 

i mental  (/;A)sLl,.k  values  in  the  1 0  to  20  W/K  range  arc  dearly 
higher  than  those  expected  from  natural  convection  alone. 
Welty  et  al.  [6]  quote  approximate  h  ranges  from  5  to  50  W / 
(nr  K)  for  natural  convection  of  air  and  25  to  250  W/ 
(m2  K)  for  forced  convection  of  air.  For  further  comparison, 
Fuller  and  Newman  [7]  and  Nguyen  and  White  [8]  proposed 
h  values  for  fuel  cell  heat  losses  to  air  in  the  range  of  25  to 
26  W/(m2  K)  for  gases  in  forced  convection  which  corre¬ 
spond  to  (fiA)^  values  in  the  range  of  7  to  10  W/K,  again 
lower  than  our  experimental  values. 

Possible  explanations  for  the  apparentdisagreements  could 
be  the  presence  of  conduction  losses  from  the  stack  to  the 
bench  and  from  the  stack  up  the  heavy  copper  power  cables, 
plus  forced  circulation  of  laboratory  air  which  results  in  heat 
transfer  further  up  the  forced  convection  scale  than  assumed 
in  the  publications  cited  above  [7,8).  The  mean  (hA)  value 
of  17  W/K  for  our  particular  stack  and  experimental  conlig- 
uration  is  considered  quite  acceptable.  This  value  may  then 
be  used  in  Eq.  ( 17)  to  predict  the  heat  loss  from  the  stack  at 
various  stack  temperatures. 

Similarly,  the  increase  in  sensible  heal  terms  can  also  be 
linked  to  heat  transfer  coefficients,  (/?A),  from  the  stack  to 
the  fluid,  /,  where  i  =  anode,  cathode,  or  water  stream,  using: 


(M)a  = 


-  o.5(  r,,.„  +  r.,.llut) 
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irochemical  model  with  an  nnsteady-siatc  thermal  model  to 
form  an  overall  transient  model.  Initially  an  overall  steady- 
state  model  was  developed  to  describe  the  steady-state  equi¬ 
librium  conditions  of  the  fuel  cell  and  then  the  model  evolved 
into  an  overall  transient  model  describing  fuel  cell  character¬ 
istics  as  a  function  of  time  until  a  new  steady-state  is  reached. 

An  unsteady-state  energy  balance  is  the  extension  of  the 
steady-state  thermal  model  by  the  addition  of  an  accumula¬ 
tion  term,  such  that: 

■Y^/d  7~S|>--t:\  _  .  .  _  . 

Ay  Cl  ^  |  f/i  [■[',.  ..  -I  i\  ^iu\s  (20) 

where  M  is  the  total  mass  of  the  fuel  cell  stack,  C  is  the 
average  specific  heat  of  the  stack,  and  dr„;,,k/dr  is  the  tem¬ 
perature  change  with  respect  to  time. 

The  estimated  value  of  MC  in  the  accumulation  term  may 
be  determined  by  using  the  summed  A/,C,  values  of  the  indi¬ 
vidual  components  of  the  fuel  cell  stack  such  as  the  graphite 
and  stainless  steel.  The  masses  of  the  individual  components 
were  calculated  based  on  the  average  density  of  the  material, 
and  the  specific  heats  were  estimated  from  data  in  the  litera¬ 
ture.  A  value  for  the  sum  of  A^C,  was  estimated  as  35  kJ/K. 
The  MC  value  may  also  be  measured  using  an  experimental 
approach  by  rearranging  Eq.  (20); 


(M 

where  the  subscript  i  represents  values  immediately  after  a 
step  change  in  power  and  the  subscript  f  represents  the  ‘final’ 
values  following  a  step  change  in  power.  An  average  value 
for  MC  was  calculated  as  approximately  35  kJ/K  which 
agrees  with  the  previous  MC  calculation. 

The  transient  model  may  be  obtained  by  coupling  the 
steady-state  electrochemical  model  and  the  thermal  model  as 
a  function  of  time  with  the  new  accumulation  term  from  Eq. 
(20)  giving: 


d  T 


stuck 


d/  MC  ^  ^thco  *7in>s ) 


(22) 


The  he;)t  transfer  parameters,  (/iA)„  as  shown  in  Table  2 
were  approximated  based  on  preliminary  data  obtained  from 
the  Ballard  Mark  V  stack.  Note  that  the  heal  transfer  param¬ 
eter  for  the  stack  is  also  shown  in  the  Tabic  for  completeness. 

Table  3  shows  a  typical  application  of  the  steady-slate  ther¬ 
mal  model  at  an  operating  current  of  20  A.  The  last  column 
shows  the  actual  experimental  data  at  the  same  operating 
conditions.  Generally,  the  model  predicts  the  experimental 
data  very  well. 


4.  Transient  model 

A  model  predicting  transient  responses  in  fuel  cell  opera¬ 
tion  was  developed  based  on  coupling  the  steady-state  clec- 


with  all  the  lerms  on  the  right  side  of  Lhe  equation  defined 
and  quantified  earlier.  This  expression  can  then  be  used  as  a 
basis  of  a  finite-element  calculation  using: 


with  a  At  value  appropriate  to  the  stack,  approximately  one 
minute  or  less  for  this  study.  The  steady-state  electrochemical 
model  imports  a  current  value  of  Tslllck  along  with  the  other 
required  system  parameters  and  then  provides  a  Kell  for  the 
calculation  of  the  updated  value  of  qc]a..  The  thermal  model 
calculates  the  new  outlet  stream  temperatures  based  on  the 
current  value  ol  so  that  lher/w,„,  term  may  be  calculated. 
The  remaining  energy  terms,  and,  </i„ss  arc  dependent 
on  TsllKk.  A  new  value  for  dT^/df  is  then  calculated  using 
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Table  3 


Typical  application  of  the  steady-slate  thermal  model  when  a  load  of  20  A  is  applied  to  the  fuel  cell  system 


Mode!  input  data 

Model  output  data 

.Stream 

Parameter 

Value 

Stream 

Parameter 

Predicted 

Experimental 

value 

value 

Intel  anode  gas 

mu 

0.0078  mol/s 

Outlet  anode  gas 

Wh-.ii.ui 

0.0041  mol/s 

0.004!  mol/s 

^  >i.in 

23.5  °C 
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25.3°C 

35  psig 

Inlet  cathode  gas 

0.004  mol/s 
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0.002  mol/s 
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23.5  °C 
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35  psig 

Inlet  water 
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23  9  °C 
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Fig.  S.  Transient  plot  for  the  start-up  of  the  Ballard  Mark  V  stack  from  0  in 
20  A. 


Fig.  7.  Transient  plot  for  the  shut-down  of  the  Ballard  Mark  V  stack  front 
60  to  0  A. 

5.  Conclusions 
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Ftg.  6.  Transient  plot  for  the  ‘toad  step-up'  of  the  Ballard  Mark  V  stack  from 
20  to  60  A. 

Eq.  (22)  so  that  a  new  7  slack  for  the  next  time  step  may  be 
calculated  using  Eq.  (23).  As  a  result,  the  transient  model 
calculates  T.uck  and  the  outlet  stream  temperatures,  Va.|1,  and 
ultimately  the  power  P  as  a  function  of  time. 

Ftg.  3  shows  the  application  of  the  transient  model  for  the 
start-up  of  the  fuel  cell  stack  for  an  operating  current  from  0 
to  20  A.  Fig.  6  shows  the  power  step-up  from  20  to  60  A  and 
Fig.  7  shows  the  shut-down  of  the  fuel  cell  stock  from  60  to 
0  A.  Each  of  the  three  plots  also  show  some  experimental 
data  for  reference.  In  general,  the  transient  model  predicts  the 
experimental  data  very  well. 
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The  transient  model  is  a  useful  tool  in  predicting  the  per¬ 
formance  of  the  PEM  fuel  cell  stack  in  our  experimental 
configuration.  For  example,  given  a  set  of  gas  feed  conditions, 
pressures,  heat  transfer  coefficients,  and  operating  current, 
the  model  will  lead  to  a  unique  solution  for  stack  voltage, 
stack  power,  and  stack  temperature.  If  the  stack  undergoes  a 
perturbation,  such  as  the  initi  al  start-up,  a  large  step  in  current, 
or  a  shut-down,  the  transient  model  will  be  able  to  generate 
transient  information  such  as  stack  temperature,  cell  voltage, 
and  power  as  a  function  of  time. 

The  steady-state  electrochemical  model  parameters,  the 
heat  capacity,  and  the  heat  transfer  coefficients  are  obviously 
specific  to  the  stack  used  and  to  the  immediate  surroundings 
of  the  stack.  These  values  would  have  to  be  recalculated  for 
another  fuel  cell  stack  system.  The  general  methodology, 
however,  should  be  applicable  to  any  PEM  stack  system. 


6.  List  of  symbols 

a  constant  in  the  Nussell  equation 

A  surface  area,  nr 

c  concentration,  mol/ cm1 
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C  average  heat  capacity,  J/ (kg  K) 

E  thermodynamic  potential,  V 

£°  standard  electrochemical  energy,  V 

F  Faraday’s  constant,  F= 96  487  C/cquivalent 
Gr  Grashof  number 

h  convective  heat  transfer  coefficient,  W/  (m2  K) 

U  change  in  enthalpy,  J/mol 

i  operating  current,  A 

I  current  density,  A  / ft2  or  ASF 

k  thermal  conductivity,  W/ (s  m  K) 

M  mass  of  the  fuel  cell,  kg 

m  constant  in  the  Nusselt  equation 

N  molar  flowrate,  mol/s 

n  number  of  cells  in  the  slack 

P  power  output,  W 

PEM  proton  exchange  membrane 

Pr  Prandtl  number 

q  energy,  W 

R  ideal  gas  constant,  R = 8.3 145  J /  ( mol 

K)  =0.08206  I  atm /(mo!  K) 

R  resistance,  H 

T  temperature,  K 

V  output  voltage,  V 

Greek  letters 

A H  heat  of  reaction,  J/mol 

rj  overvoltage,  V 

f  parametric  coefficient  in  the  fuel  cell  model 

Subscripts 

a  anode 

act  activation 

c  cathode 

cell  proton  exchange  membrane  fuel  cell 

cons  consumed 

elec  electrical 

g  gas 

hum  humidification 

in  in 


int  internal 

1  liquid 

loss  loss 

ohmic  ohmic 

out  out 

prod  product 

rxn  reaction 

stack  fuel  cell  stack 

sens  sensible 

theo  theoretical 

w  water 

Superscripts 

sat  saturation 

*  at  the  catalyst  interface 
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